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THE  FLOW  FIELD  SURRODIOINQ  A MOVINQ  VESSEL 


Introdttctioa 

The  calculation  of  a flow  pattern  is  the  first  step  in  ever/ 

theoretical  explanation  of  pressure  fields.  Theories  of  pressure 

fields  are  nunerous;  In  a sense,  therefore,  the  substance  of  the  present 

repo^  is  not  new.  Yet,  In  aost  of  the  work  that  has  come  to  the 

author *s  attenti<m,  the  emphasis  Is  so  preponderantly  on  matters  other 

than  flow,  and  discussion  of  the  water  motions  Is  so  sketchy,  that  more 

wox^  with  slxigolar  stress  on  flow  seemed  greatly  needed.  Furthermore, 

(a)(*) 

in  one  of  the  best  reports  on  pressure  fields  , which  use&the  model  of 
a line  source  and  a line  slnk^  transverse  velocities  of  the  water  are 
omitted  from  consideration.  Extension  of  the  work  to  fill  that  gap  Is 
required. 

Hence,  the  purpose  of  the  present  account  Is  twofold:  (1)  It  Is 

to  review  the  hydrodynamics  of  the  problem,  outllrlng  the  methods  that 
have  been  found  fruitful  and  worth  knowlaig.  In  attempting  this  we  shall 
not  be  bound  by  the  usual  inhibitions  of  the  writers  of  technical  reports 
and  present  all  the  material  necessary  for  a coherent  picture->-even 
though  most  of  it  Is  available  in  the  literature.  For,  while  the  other 
procedure  saves  the  writer’s  time.  It  often  wastes  the  readers’.  (2) 

We  desire  to  provide  the  numerical  material  concerning  the  theory  of  flow 
that  seems  to  be  still  missing. 

Ttie  disturbance  created  by  a moving  ship  can  be  simulated  by  a 
variety  of  mathematical  models,  best  peiiiaps  by  combinations  of  sources 
{^)  Small  leHers'  above  thie  line  refer  to  kotes,  which  are  on  page  5^. 
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(b) 

and  sinks.  T«ro  such  models  are  here  investigated.  The  first  consists  of 
a point  source  and  a point  sink  placed  slightly  below  the  svirfaee  of  the 
" water.  The  second  model  is  c(»aposed  of  a line  source  of  constant  strength 

/ near  the  bow  of  the  ship  and  a similar  line  sink  near  Its  stem.  A point 

source  discharges  water  in  radially  symnetrio  fashion  and  at  a uniform 
rate  about  a point,  a point  sink  absorbs  It  in  similar  fashion.  Henee- 
forth,  where  it  will  be  unambiguous,  we  will  say  "souroes"  to  mean  "sources 
and  sinks",  etc.  The  line  source  may  be  thouj^t  of  as  a long  pipe 
uniformly  perforated  on  its  lateral  surface  and  exuding  water  at  a 
Wiifcrm  rate.  A point  souree  is  said  to  have  strength  m if  it  emits  a 
volume  of  water  equal  to  4'rrm  ft^  per  second,  i.e.  if  it  causes  a flow  of 
» ft3  per  second  into  every  solid  angle  of  unit  magnitude.  The  strength 
* of  a line  souree,  f(f),  is  la  general  a function  of  the  coordinate 
of  a point  in  the  source.  We  shall  define  it  in  such  a way  that: 

jf(|  )<J|  ■ n,  (1) 

the  total  strength.  In  this  ease,  4 ir  m is  again  the  total  volume 

of  water  emitted  per  second.  In  our  application  f ( 'I ) is  constant  and 

(c) 

equal  to  m/X  provided  JL  is  the  length  of  the  line  source.' 

To  adjuet  OuT  mOuvIb  tO  actuol  COuditlOuo,  vuo  da  Vice  cf  linagea 
is  employed.  These  will  be  so  disposed  that  the  flow  is  parallel  to 
the  surface  and  to  the  bottom  of  the  water,  as  explained  below. 

^ Vlseosity  effects  are  not  ocmsidered. 

Mathematical  Preliminaries. 

Since  the  flow  outside  the  sources  and  sinks  is  everywhere  irrota- 
tlonal,  velocities  can  be  derived  fr<»n  a velocity  potential  ^ with  the 
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For  a point  source 
For  a line  souroe 


m 
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(2) 

(5) 
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'V'  S V ^ 

'X/' 

4> 

4’  = iUUiS. 

r being  the  distance  froa  the  source  element  at  ^ to  the  point 
where  ^ is  required. 

In  addition  to  the  velocity  potential  we  shall  need  the  stream 
function,  . The  ordinary  stiream  function  is  useful  in  connection 
with  problems  of  two*dimensional  flew.  It  is  defined  as  'le  rate  of 
flow,  from  right  to  left,  across  the  curve  OP  (see  Fig.l).  In  the 

absenee  of  sources  and  sinks,  this  rate  of 
flow  is  Independent  of  the  curve  and  depends 
only  on  the  coordinates  of  the  point  P.  l^us, 

^ -V"U.7). 

The  water  displaced  by  a moving  ship 
does  not  flow  in  two  dimensions,  even  approx- 
imately. A more  reasonable  first  apprcxlma-  Figure  1 

tion  would  consider  axisyms^tric  motion,  such  aa  results  from  translation 
through  the  water  of  a cigar-shaped  body  along  its  axis.  A modified 
stream  function,  designed  to  describe  axi«yB5«et-riQ  motions,  was  Intrc-^ 
duced  by  Stokes.  Let  the  axis  of  symmetry  be  X and  use  the  set  of 
cylindrical  coordinates  presented  In  Fig.  2.  Select  a point  A on  the 
axis,  conneet  it  by  a curve 
to  P,  the  point  where  the 
stream  function  is  to  be 
evaluated.  Finally,  con- 
struct a stirfaoe  of 
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revolution  by  rotating  that  curve  about  X.  The  flow  to  the  left 
across  this  surface.  In  terms  of  Stokes*  stream  function  ^ » is  2.1T^ , 
Prom  this  definition  it  is  clear  that  ^ will  In  general  depend  on  the 
coordinates  of  P and  A,  as  well  as  on  the  eux*ve  connecting  these  points. 

In  the  absence  of  sources  and  sinks,  however,  the  curve  is  irrelevant 
and  If,  furthermore,  there  is  no  flow  across  the  X-axis, 

the  fluid  crossing  PA  will  also  cross  FB,  and  the  position  of  the  point 
A does  not  affect  provided  A is  on  the  axis.  Both  of  these  conditions 
are  satisfied  in  the  problem  under  study.  Hence  we  conclude  that  V" 
is  a function  of  the  coordinates  of  P alone.  In  the  rraalnder  of  this 
report  denotes  Stokes*  stream  function.  It  follows  from  the  definition 

Just  given  that  “dp  ' ^ 9x^ 

Note  that  the  physical  dimensions  of  Stokes*  stream  function  are 
idiereas  those  of  the  ordinary  are  . Since  there  is  no  flow 

across  a curve  of  constant  the  equation  constant  defines  a 
streamline.  Herein  lies  the  principal  lnqportance  and  the  usefulness  of 
the  stream  function. 

The  fluid  emitted  by  a point  source  at  o through  the  spherical 


cap  ?P  {of  Pig.  3)»  fres  right  to  left,  is  - mill 
ilL  being  the  solid  angle  subtended  by  the  cap 
at  0.  But  jfi  • 2 7T’(l-cos9).  Hence,  by  defin- 
ition, va(  coiS  ^ l). 

It  the  source  is  located  not  at  0 but  at  some 
point  x^  on  the  X-axis, 

„ _ l\ 
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We  shall  now  find  the  stream  function  generated  by  a line  source 
extending  along  Consulting  Fig.  4 we  find  that 


Here  corre8p<»d8  to  the  left,  03  to  the  right  end  point  of  the  line 
source.  If  we  denote  the  distances  of  these  end  points  from  P by 
and  T2  we  have  finally 


(5) 

There  remains  the  simple  task  of  determining  the  stream  function 
for  uniform  flow  of  an  infinite  body  of  fluid.  If  its  velocity  to  the 
rli^t  is  7,  then  clearly  the  flow  to  the  left  through  a circle  of  radius 
^ perpendicular  to  X,  idiieh  by  definition  is  2ir  , equals  -ir^V. 
Therefore  P^V. 

The  flCw  pattern  produceu  by  a vingie  souz*oe  is  a radial  outward 
current.  A source  and  a sink  together  generate  streamlines  identical 
with  the  magnetic  force  lines  surrounding  a bar  magnet  with  two  poles. 

* HSither  model  bears  any  relation  to  the  flow  pattern  around  a ship. 

But  if  a source  and  a sink  are  embedded  in  a stream  of  constant  velocity 
V,  the  resulting  flow  simulates  the  situation  in  whloh  water  flows 
parallel  to  the  keel  of  a stationary  boat.  There  will  be  two  points 
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where  the  velocity  imparted  to  the  water  by  the  source  sink  system  Is 
zedueed  to  zero  by  the  velocity  of  the  stream:  These  are  stagnation 

points  and  correspond  to  the  bow  and  the  stem  of  the  boat. 

It  Is  Indicated*  therefore*  that  a useful  for  present  purposes 
is  a superposition  of  a stream  function  of  seme  souree*slnk  model  with 
one  representing  unifona  flow:  Vs  ” ^ p K (6) 

For  w»jhave  chosen  different  forms*  and  these  will  be  discussed 
In  due  course. 

l!he  models  described  thus  far  simulate  the  motion  of  cigar  or  egg* 
shaped  objects  throued^  an  Infinite  medium  In  the  direction  of  their 
axis  of  (highest)  symmetry.  They  ignore  the  presence  of  boundaries* 
which  destroy  axial  synnetry.  One  of  the  effects  of  surface  and  bottom 
Is  to  enforce  a flow  which  Is  wholly  parallel  to  these  surfaces  in  their 
Immedlavd  selghbozhood.  A formal  method  for  guaranteeing  this  in  our 
theory  is  available  through  the  use  of  images*  as  follows: 

Assume  tiie  water  to  have  a depth  h»  Theory  takes  the  water  to  be 
infinite  In  all  directions*  but  singles  out  a hypothetical  "slab"  of 
thickness  h*  the  horizontal  surfaces  of  which  correspond  to  the  real 
fliirfafia  aisd  the  real  bettes.  A scares  anywhor-e  within  the  "siab"  will 
cause  a vertical  flow  across  the  surface*  since  the  source  knows  no 
boundary  conditions.  But  if  another  source  of  equal  magnitude  is  placed 
at  the  image  point  of  the  first*  beyond  the  surface*  Its  vertical  flow 
will  cancel  that  of  the  primary  source.  Both  of  these  sources*  however* 
send  water  across  the  aea  floor.  Ti'  raiedy  this  fault*  both  must  be 
neutralized  by  a eonstruotion  of  images  beyond  the  bottom  surface. 
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These  must  then  be  reflected  again  In  the  top  surface,  and  so  on,  until 
an  infinite  set  of  sources  in  a vertical  line  is  produced,  which  serves 
to  satisfy  the  boundary  conditions  with  respect  to  vertical  flow^ 

If  the  primary  source  lies  a distance  'f  below  the  sux*faoe  and 
we  take  the  origin  within  the  surface,  images  must  be  constructed  at 

In  particular,  if  the  prlxoary  source  is  situated  in  the  surface  Itself 
images  occur  at  - 2 Jh. 

The  presence  of  images  maJtes  the  treatment  of  certain  aspects  of 
the  problem  difficult o This  is  especially  true  with  respect  to  the 
details  of  flow  near  the  hull  of  the  vessel;  for  instance,  tc  find  the 
exact  streamline  which  defines  the  boat  contour  in  the  presence  of  all 
images  is  a tedious  task.  Fortunately,  it  need  not  be  undertaken,  for 
the  flow  near  the  vessel  is  controlled  pred<MBlnantly  by  the  primary 
source  (or  sources).  To  obtain  it  we  Ignore  the  images  and  treat  the 
problem  as  though  we  had  axial  sysBaetry.  But  the  computation  of  the 
water  motions  near  the  bottom,  which  are  the  result  of  the  source  and 
many  lomigeB,  must  take  account  of  the  Infinite  array  of  images. 

we  now  present  the  special  features  of  the  two  models  chosen  for 
careful  study. 
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Point  Source  Model 

Our  method  for  both  models  is  to  conqpute  the  velocity  on  the  sea 
floor  caused  by  a stationary  source  distribution  lying  in  water  which 
is  stagnant  at  great  distances.  The  ship  speed  enters  the  computations 
only  thru  the  strengths  of  the  sources  which  are  proportional  to  V. 

This  steady  velocity  is  tdie  same  as  that  which  would  be  observed  Instan- 
taneously when  the  distribution  moved  by  with  the  speed  of  the  ship, 

provided  of  course  that  the  vector  distance  between  the  distribution 

(d) 

and  the  point  of  observation  were  the  same  in  both  eases. 

The  present  model  approzlmates  the  ship  by  a point  source  near  the 

bow  and  a point  sink  near  the  stem.  The  water  surface  and  the  sea 

floor  are  ai^raximated  by  rigid  horizontal  planes.  The  boundary  con- 

dlti<ms  - that  there  be  no  flux  across  these  planes  - are  satisfied  by 

constructing  image  sources  and  sinks  as  described  on  pages  7 &nid  8. 

Let  the  water  surface  be  the  X7  plane;  the  sea  floor,  the  plane  ^--h  , 

and  let  the  source  and  sink  be  a distance  f below  the  surface  and  a 

(e)  ” 

distance  2JL  apart  along  the  z axis. 

First  we  calculate  , caused  only  by  the  source  and  its  images, 

fc-t-  * _ « * /*  I ■ 'n  \ — / rs  _ J.  \ ^ 

uu  Bva  XJ.WWX-  «kXMii(s  jcecx  x.xjt*o,  ~ ^ , xx-wui 


xeO  to  500  feet.  To  do  this  we  must  sum  the  individual  contributions 
from  the  source  and  Its  Images.  Ohe  general  z-eoordinates  of  the 
images  have  been  given  in  Eqs.  (7)t 

2yL  4-  = + = W-,  - g 

- 0 \ I ^ ^ ^ I ^ • * 

-(Z  i U 4.  f) 


2 

- 


is  1,1,1'"®^ 
s I ) 2.  ' • ' 
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Prom  Eqs.  (2),  (3),  and  (7) 


~v^  = yy7^y-a,) 


oo 

^ ("if  ^ 


(9) 


wnero  y - {X'^c^j  t*  ( .-  ( c - v-^  j ^ wiwrie  u,  e^  ) 

are  the  coordinates  of  the  source  and  its  images;  s«  1,2, 3, 4-. 


Also, 


• -^  = ^£f 

I.C  I 


(rrr 


H 


We  must  nov  ecMBpute  the  sum. 

under  certain  conditions  the  sum  of  a series  may  be  approximated 
(and  its  convergence  simultaneously  proved)  by  Integrals.  If  we  denote 
by  u^^  the  general  term  of  the  series  and  by  u^O  the  function  formed 
by  allowing  1 to  take  all  values,  we  have 

so 

^ W.^  O i i- 


u.«>  + 


or 

<o 


t/L 


f^ichsver  limits  the  sum  sore  nai-rcwiy.  The  conditions  on  these 
inequalities  are:  (1)  The  Integrals  must  exist;  and  (2)  beyond 

i,e  'S  ^ u.(D  isust  be  monotonically  decreasing,  . This  can  be  seen  from 

0O 

Pig.  6.  The  area  represented  by  the  integral  fu.(L)6i  on  the  left 
side  of  the  Inequality  is  indicated  by  shading  thus: 
the  right  by:  , tui^i 


Figure  6 
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Then  as  an  approximation 


if 


•h 


(10) 


One  of  the  terns  In  the  sunioatlons  constituting  (eqc  9) 


18 


CO 


vy\  (X- 


where  p**  h 

80  the  Integral  we  need  In  Eqs.  (10)  Is 


o ! ^ +.  C0«s1"^v»f> 

where  ^ = a - . For  seciuences  Ho,  1 and  2,  d'5=  ~2Vi«Ju  ; for 

sequctncea  Nos.  3 and  . For  1 and  2 the  limits  become 

— os*  iKv\6  . For  3 and  4-  the  limits  bec(»Be 

oo  a.vvJ  jf- 

Since  the  point  of  evaluation  most  lie  within  the  water  all  the 
conditions  on  the  integrals  in  Eqs.  (10)  are  satisfied;  this  can  be 
seen  most  readily  from  Pig.  3.  As  ^ Increases  the  distant  images 
have  a proportionately  greater  effect  on  the  velocity^  so  that  the  series 
converge  less  rapidly. 
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Then  for: 

vn(x-o^) 

> 

S»3,¥  J,1to 


oo 

/ ^ - yx\CX-<^)  \^  t 

~ h L 

a JXiLh:r  ^lf^  i 

L 


(11) 


To  ooaoqpute  'V^  asisg  Eqs.  (11)  we  must  select  values  of  h,  m, 
and  f . To  be  able  to  coapare  our  result  with  8<sBe  of  those  In  NOL 
Report  Ho.  504»  we  have  taken  h ■ 54  ft.»  the  shallowest  depth  treated. 
According  to  Lunde«  op.  clt.  p.  38,  m * » AV , where  A is  the  transverse 

— 4TT 

mldseetlon  area  and  V Is  the  ship  speed;  and  f Is  the  depth  of  the 
centroid  of  A.  The  average  specifications  of  the  four  ships  for  irtileh 
the  experimental  pressure  signature  Is  given  In  Reference  (2)  (plates 


82,  87,  88,  93)  are: 

bean  B 

draft  D 

length  L 

tonnage 

displacement  A 


60  feet 
24.5  " 
441-7A2  feet 
6838  tons 
214,000  ft5 


If  the  transverse  mldseetlon  were  rectangular,  A would  equal  BD;  and 
If  It  were  seml-^lllptlcal,  A would  be  (Tr/4)  BD.  Its  shape  is 
probably  between  these  two:  hence  we  take  A ■ (0.9)  BD.  Then 
A s (0.9)  (60)  (24.5)  * 1320  ft2-  Take  m • 2000  ft3/sec.  For  the  ships 
under  ocmslderatlon  this  strength  corresponds  to  a speed  of  11.5  knots. 
The  velocity  Induced  by  a ship  of  different  cross-section  or  moving 
with  a different  speed,  can  be  obtained  froaa  our  results  by  one 


0OT?PIlX5m.ftI. 
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multiplication.  Again,  if  the  mldBectlon  is  rectangulsor  f ^ 12.25 
feet.  Since  it  is  somewhat  elliptical,  we  take  f « 10  feet, 

and"\'^  at  any  olher  point  on  the  sea  floor  may  be  derived 
from'v^  along  the  keel'^line  by  noting  from  Pig.  7 that  for  points 
on  a circle  of  constant 


'Vi-  = VC  JL 


> 


- v^Jd. 


(12) 


Y 


With  the  use  of  these  relations,  "v^  and^  were  obtained  as  functions 
of  X for  y ■ 0,  5,  15,  25,  40,  50,  75»  100  feet.  All  these  data  are 
presented  in  Figs.  15  and  l6. 

The  quantity  X.  was  determined  by  requiring  that  the  points  on  the 
sea  floor  beneath  the  keel  at  which s o be  the  same  distance  iqjart 
for  both  the  point  and  the  line  modeX,  This  gave  s itv  feet. 

Having  found  jt  and  knowing  by  inspection  the  signs  of  the  velocities 
contributed  by  the  source  and  the  sink,  the  total and*\xj  along  each 
line  t|  « constant  was  obtained  from  the  graphs  by  appropriate  additions 
and  subtractions.  The  contours  of  constant and'^vzj  on  the  sea  floor 
beneath  the  ship  are  shown  in  Figs.  18  and  19. 
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To  get  a roizgh  estimate  of  the  size  of  the  eloeed  stresm-surface 
the  dimensions  of  a simple  souree-slnk  Fankine  solid  were  computed  for 
jrvi  ^ ifooo  (to  include  the  fix*st  surface  images).  The 

appropriate  theory  is  given  in  Hilne->Th<HBoon,  pp  411-412.  F:  <»b  this 
calculation:  B * 57*  D " 28.5,  L ■ 582  feet.  Since  the  solid  is  blunt 
and  the  ship  sharp,  making  the  equivalent  solid  shoi*ter,  eosq;>an8on 
with  the  average  dimensions  of  the  xeal  ships  treated  (p  5)  shows  that 
the  m and  1 we  have  used  are  substantially  correct. 

The  streamlines  and  stagnation  points  la  the  XY  plane  near  the 

sources  are  shown  In  Fig.  20,  and  the  l<mgitudlnal  and  transverse 

sections  of  the  closed  streamsurfaces  of  our  source  distribution  are 

(f) 

sketched  in  Fig.  21.  Bote  that  the  solid  is  divided  into  two  separate 
parts,  sysanetrical  about  the  XY  plane.  The  length  of  these  solids  is 
about  374  feet;  the  beam,  less  than  57  feet;  the  draft,  greater  than 
28.5  feet.  Thus,  the  actual  solid  does  not  approximate  the  dinensiinut 
or  the  shape  of  a ship  as  well  as  the  simple  source-sink  solid.  If  one 
were  able  to  fit  a streamsurface  exactly  to  a ship  hull  the  velocities 
derived  should  be  exact  within  t^e  limitations  of  the  other  approxlaa- 
tlons  made.  But  since  the  values  we  find  for  the  velocities  do  not 
differ  appreciably  from  those  obtained  with  the  use  of  a line  source 
model  whose  streamsurface  fit  a ship  unusually  well,  we  may  conclude 
that  this  "goodness-of-fit”  eriterionp  based  upon  an  examination  of  the 
dividing  streaiallns,  does  not  provide  on  n priori  test  of  the  validity 
of  velocity  calculations. 
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Since  the  results  of  the  present  section^  vrhich  are  based  upon 
the  point  source  model,  are  not  very  different  from,  and  indeed  tend 
to  confirm,  the  results  of  the  next  section,  no  separate  discussion 
of  the  present  data  will  be  given.  It  is  Interesting  to  note,  however, 
that  the  maxima  and  minima  are  more  closely  spaced  for  the  point  model, 
and  that  the  velocity  gradients  are  somewhat  steeper  near  these  maxima 
and  minima  than  In  the  line  model.  The  present  model  was  used  to 
obtain  evidence  about  the  trustworthiness  of  the  different  ways  of 
simulating  ships,  and  it  is  gratifying  to  note  that  computed  velocities 
are  so  remarkably  independent  of  the  choice  of  model. 

Linear  Source  Model 

A line  sink  and  a line  source  of  equal  strength  m and  of  equal 
length  4.  are  placed  along  the  X-axis,  as  shown  in  pig.  8. 


The  X-axis  Is  In  the  surface,  the  p-axls  extends  at  right  angles  to 
It.  Our  analysis  here  follows,  and  In  many  respects  repeats  the  work 
of  NOL  Report  No. $04.  The  stream  function  for  this  system  is  given  by 
Eq*  (3)  which  must  now  be  applied  to  each  of  the  two  line  sources,  the 
first  with  negative  m.  Hence, 
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To  simulate  fluid  motion  past  a boat^  we  superimpose  the  constant 
velocity  • V*  (5)  then  leads  to  the  expression  for  the  stream 


lines: 


-ir  X p^V  a constant 


(15) 


In  this  approximation  images  have  been  ignored  and  boundary  conditions 
are  not  satisfied.  It  is  fairly  certain,  however,  that  the  lower  lialf 
of  the  axially  symoetrie  streamline  pattern  represented  by  Eq.  (15) 
will  portray  reasonably  well  the  flow  in  the  neighbor nood  of  the 
primary  source  even  though  the  flow  is  perturbed  by  the  images* 

Eq.  (15)  nay,  therefore,  be  used  to  detemine  the  streamline  which 
contains  the  two  stagnation  points.  Often  it  is  called  the  "dividing 
rtreamllne"  because  of  its  branching  at  the  stagnation  points.  By 
proper  choice  of  parameters,  this  dividing  streamline  can  be  made  to 
fit  a boat  very  well  (see  Pig.  22),  so  we  ragy  expect  our  velocity 
results  to  be  satisfactory. 

Par  ahead  of,  and  far  behind  the  boat,  this  streamline  must  clearly 
lie  on  the  X-axis.  That  condition  is  xaet  if  in  Eq.  (15)  we  put  the 
constant  equal  to  zero,  for  the  left  side  is  identically  .?ero  at 


^ A 4 Mm  m 4 A 

M.L  vr  U Aaaw 


must  be  given  everyidiere  by  the  condition 

^ + -V  ?*V'  = 0 (16) 

To  repzvseat  a real  boat,  we  adopt  the  recommendation  of  NOL  Report  #504 
and  take  jl»tZS  ^ =8  3/5“  •fi,*' 

Eq.  (l6r^a  plotted  for  these  values  in  Fig.  22^  The  fit  of  the 
present  vodel  Is  exeellent.  2t  eorxespoads  to  a vessel  4yo  ft.  long. 
Perhaps  it  Is  also  of  interest  to  recognise  that  cur  choice  of  constants 
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implies  a source  strez^gth  m of  5360  ft  '/sec  for  a vessel  moving 
with  a speed  of  10  knots,  that  is,  an  exchange  of  about  67,500  cubic 
feet  o)?  water  per  second  between  source  and  sink.  We  now  proceed  to 
the  calculation  of  flew  velocities. 


A.  Velocity  along  the  direction  of  the  ship»s  motion. 

Since  the  lodges  remove  the  axial  symmetry  of  the  problem,  the 
advantage  of  cylindrical  coordinates  disappears  and  it  now  becomes 
preferable  to  use  Cartesian  coordinates.  Ihe  Z*axls  is  taken  vertical » 
and  X is  in  the  direction  of  V within  the  water  surface,  ihe  present 
task  Is  to  find  v^^ 

According  to  Eq.  (4)  . ^ 

, p.  V(x-tJ’+y+z‘ 

r 

and  our  model  Is  specified  by 


, and  a partial  integration  leads  to 

- (O’J'*? 


In  the  present  case  f (f  ) vanishes  at  the  limits  of  integration. 

Purthe^ox,.  ^ f (y,X.) * f (f,  X3)  - / (f , V)3 


being  Dirac's  singular  function.  Consequently 

jn_  (J, ! !_  ^ 

. ’’i  ''m 


■vifp)* 


l Vr 


!r) 

"4. 
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where  the  r*»  are  the  (positive)  distances  of  the  point  P from  the 
foor  ends  of  the  two  line  sources.  Interestingly,  due  to  our  line- 
source  distribution  la  identical  with  the  potential  resulting  from 
four  point  sources  of  magnitude  f*  situated  at  the  four  end 

points  of  our  lines. 

At  tills  stage  all  Images  are  Introduced,  and  v^  will  be  computed 
for  a point  P on  the  sea  floor.  The  situation  is  depicted  in  Pig.  (9). 
It  is  necessary  to  consider  the  effects  of  four  infinite  vertical 
sequences  of  point  souz^es,  occurring  in  the  figure  below  the  numerals 
1,  2,  3 and  4,  Thus  (rj^^  la  defined  in  the  figure) 
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The  result  as  stated  is  valid  for  a 


•a 


= 1,  2,  3.  4. 
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These  sums  converge  too  slowly  for  direct  evaluation.  However,  each 

can  be  converted  into  another  rapidly  converging  series  by  a device 
(g) 

of  Made  lung  whose  interest  was  in  an  electrical  analog  of  the 

present  problem.  He  has  shown  that  , 

Here  d is  the  horizontal  distance  of  P from  the  vertical  row  of 

points  under  consideration:  d;  - 

(h)  J 

H is  the  Hankel  function 

point  F whose  vertical  position  is  midway  between  two  sources,  and  not 

in  the  line  of  sources.  Constant  terms  which  eaxicel  when  Eq.  (I?)  is 

formed  have  been  omitted.  In  a line  of  sources,  l.e.  when  dj  = 0,  L 

takes  on  the  limiting  value  L(o)*»  I.963. 

When  the  point  p lies  in  the  surface,  we  are  also  led  to  Eq.  (17)  , 

but  with  an  altered  meaning  of  the  The  summations  then  yield  a 

result  similar  to  Eq.  (l8)s  however,  all  terms  in  %.  are  now  positive. 

n 

This  modified  function  will  be  called  L^;  it  differs  from  L only  for 
small  values  of  the  ax^ument.  Both  L and  have  bean  computed  and 
graphs  are  shovnn  in  Fig.  23*  The  function  Lo  bee<xae8  00  for  zero  argument. 
In  conclusion,  then. 


X = 


m 

IT 


As  it  stands,  this  formula  is  true  for  a point  on  the  bottom;  in  the 
surface  the  same  foxv;ula  holds  but  L must  be  replaced  by  Lo  evez^where. 
With  the  data  adopted  in  the  preceding  section,  m/ajt  * 2.35  ft  V/h. 
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The  ratio 


•vr 


is  therefore  directly  given  by  Eq.  (19).  This  ratio 
•was  computed  in  NOL  Report  #504  and  plotted  in  the  form  of  contour 


graphs,,  which  we  have  checked.  As  supplementary  material  we  present 
here  some  calculations  of  the  variation  of  in  plane  vertical 
sections  p\arallel  to  X. 

Fig.  24  shows  this  ratio  in  the  keel  section  of  tlie  boat,  where 
it  takes  its  largest  values.  The  line  source  and  the  schematic  form 
of  the  boat  are  sketched  in.  At  the  bottom  of  a 40  ft.  chaxmel  *^0 
•vessel  is  preceded  by  a positive  current,  l.e.  a flow  in  tlie  direction 
of  the  vessel,  which  is  2^  of  V at  a distance  of  200  ft.  in  front  of 
the  bow.  As  the  bows  passes,  this  rises  to  8$^  and  then  falls  rapidly 
to  negatl've  values.  Under  the  left  end  of  the  line  source,  when  about 
lA  of  the  ship*s  length  !ias  passed,  ■ reaches  its  maxlmam  negative 
•value  of  nearly  0.14,  corresponding  to  2.3  ft/see  for  a ship  moving  at 
10  knots.  This  flow  diminishes  slightly  but  continues  at  more  than  the 
ratio  0.1  until  about  3A  of  the  ship's  length  has  passed,  when  the  flow 
returns  to  its  positive  sign,  the  water  tlllijog  in  behind  the  vessel. 

The  pattern  drawn  i:i  Fig.  24  is  synaoetrio;  it  repeats  itself  at  negative 


values  of  Tt<^e  rise  of 


near  x » 0 appearing  on  this  diagram 


occurs  also  In  the  point  source  model  treated  in  previous  section.  It 

-vr_ 

is  absent  at  greater  depths  as  is  shown  in  graph  25  which  gives  “ 
at  6o  ft.  Here  •the  maxlTBum  value  of  the  ratio  is  .065;  this  indicates 
that,  below  the  keel,  flow  velocities  decrease  store  rapidly  than  ^ 
but  not  quite  as  rapidly  as  lA^« 
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Figures  24  and  25  also  show  — — on  the  bottom  of  40  ft.  and  6o  ft. 

V 

channels,  a dlstanoe  of  50  ft  away  from  the  keel  section,  beam  of 

the  ship  model  used  Is  about  70  ft,  hence  the  flow  here  considered  is 
no  longer  under  the  ship.  It  is  evident  that  the  lateral  diminution 
of  X is  very  rapid,  particularly  at  the  greater  channel  depth. 

This  conclusion  can  also  be  read  from  the  contours  of  NOL  Report  i$(304. 

In  interpreting  the  last  two  graphs  It  should  be  remembered  that  -v; 
has  the  same  value  for  y and  - y. 

The  behavior  of  v^^  to  the  surface  is  not  particularly  interesting, 
although  It  llli»trates  some  features  of  our  model.  In  Fig.  26  a graph 
of  to  the  surface  Is  given.  The  ratio  Is  infinite  at  the  exact 
end  points  of  the  line  source.  Beyond  the  exterior  end  it  descends 
very  fast  to  small  values,  orosstog  the  value  1 at  a point  ;)u8t  outside 
the  line  source  (about  ^ tt,  outside  to  this  instance).  O^s  point  is 
the  bow  of  the  boat,  a stagnation  point  where  the  water  moves  with  the 
velocity  of  the  ship.  The  part  of  Fig.  26  to  che  left  of  the  stagnation 
point  has  no  physical  significance  since  it  is  excluded  by  the  body  of 
the  vessel. 

At  large  distances  ahead  of  the  boat  the  velocity  is  given  by  a 
simple  formula.  Let  the  distance,  measured  from  the  bow,  be  f and 
suppose  that  f ^ h,  and  f • One  may  then  approximate 

L,  **  Irt  ^ . 

and  (.orfae.)  « *,3^  + J 

^ fX  jE 

" IT"  iT-fTHu 
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The  dependence  on  1/ti  is  interesting 


B.  Velocity  transvei^se  to  the  dirv^ction  of  the  ship’s  motion. 
We  have 


9<j!>  _ 


kf^- 


f+R"; 


Vx 


2 2 2 2 

for  a single  source  of  strength  density  f ( ^ )•  Here  R « x 4 y t z 


Thus  our  prljaary  source  yields 


% 


(19) 


-X 


,> 

by  an  e,l«nentary  Integration.  Here  again,  the  X*azis  is  in 
surface  and  extends  in  the  direction  of  the  boat;  z is  vertical;  Xg 
is  the  coordinate  of  ah  end  point  of  the  primary  source  system;  r^  is 
the  distance  fr<»i  one  of  these  end  points  to  (z,y,z),  where  Vy  is  to  be 
evaluated. 

we  now  introduce  the  images  depleted  in  Fig.  9.  Each  mEdces  a con* 
tributicn  liks  Eq.  (19)  but  with  its  own  z(-)  and  r(^J.  Hence  we  find 


(20) 


where  only  the  st^mation  for  end  points  has  been  written,  and 

^ (in -hi) 

Ohe  summations  encountered  in  Eq.  (20)  converge  rapidly  and  can  be 
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computed  directly. 

We  define  the  fttcietion 


T 1 b)  — E b % -h/ J (21 ) 


T is  a pure  number.  In  terms  of  it 


f-lpr  LW- 


Evaluation  of  T proceeds  by  dir  ?t  suBBiation  up  to  n 

(t^O" » 


t*  such  that 
<23 ) 


^e  remainder  is  eooQmted  by  converting  the  sm  into  an  integral: 
With  the  use  of  elementary  formulas  this  becomes 


[0*^+  4 (t+O**] 


and  under  condition  (23) 

“ TotTo* 

^ # 

The  foregoing  considerations  permit  the  eomputatlcn  of  *r  wit-hout- 

mueh  labor,  we  present  the  results  in  Figs.  27  and  28.  Values  for 
intermediate  b*s  are  easily  found  by  inteirpoletlon  between  the  three 
curves  plotted. 

Eq.  (22)  shows  that  Vy  • o below  the  keel  of  the  boat  (y  « o).  In 
every  transverse  section  it  reaches  a (on  the  bottom)  some 

distance  from  the  ship  and  then  goes  gradually  to  zero.  Our  numerical  work 
shows  that  this  maximum  eomes  at  about  y ® h.  Also,  v~  differs  In  sign 
on  the  two  sides  of  the  keel  line^  as  expected.  The  quantity  In  the 
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bracket  of  Eq,  (22)  changee  sign  vfhen  x is  reflected  ir-  the  point 
midway  between  the  sources.  Hence  the  transverse  flow  pattern  under 
the  aft  half  of  the  boat  Is  opposite'  to  that  under  the  fore  half. 

Details  are  given  in  Figs.  29  and  30.  In  a channel  of  depth 
40  ftj  the  transverse  vel»>clty  reaches  a maximum  of  about  0.1  V. 

At  6o  ft.  this  drops  to  about  O.03  V.  The  maximum  always  appears  under 
the  center  of  the  sourt^c;  line.  For  y » 2h  the  flow  is  less  than  its 
Haxiimnn  by  20  to  30$^.  ^e  total  flow  Is,  of  course.  In  accord  with 
the  8ouro8-slnk  model  here  employed  and  is  represented  schematically 
In  Fig.  10. 

Figure  10 

To  eoiqplete  the  picture  we 
have  plotted  In  Fig.  31 
curves  of  constant 
on  the  bottcMa  of  a 40  ft. 
channel. 

Heasurements  of  flow  velocities  beneath  moving  vessels  have  not 
Cwiie  to  the  authors*  attention.  The  theory  oen,  therefore,  be  tested 
only  In  an  indirect*  way  by  egaimtlng  pressures  find  c<HBparlng  them  with 
observed  data.  This  is  done  in  Appendix  Xix. 
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Pi*oof  that  the  Computed  and  Observed  Velocities  art?  Equivalent . 


The  proof  of  the  theorem  stated  on  paa;e  2 la  illustrated  step  by 
step  below,  from  the  hypothetical  situation  with  which  the  theory 
deals  to  the  actual  sltuatlono  F is  a point  fixed  relative  to  the 
source  dlstrlbutloa  and  P,  the  point  at  which  the  ccmputatlons  (or 
n^asurements ) are  made»  First,  the  fictitious  case: 

Figure  11a 


velocity 
of  P ■ o 


stream  velocity 
at  CO  * o 


Call  the  velocity  of  water  at  P This  is  the  velocity  we  compute. 

Next  consider  the  source-sink  distribution  in  a stream  of  velocity-?. 

Figure  lib 


velocity 
of  P • o 


stream  velocity 
at  oo  s -V 


The  water  velocity  at  P is  now  v^  - V, 


If  the  point  P (that  is.  the  sea  floor)  also  moves  with  the  velocity 


-Vi 


velocity 
of  P • o 


velocity  of  P • 


F 

-V 


Figure  11c 


stream  velocity 
at  oo  s -V 
^ 


The  water  velocity  at  P (iwjasured  relative  to  P)  is  (v^^'V)  - V* 
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This  situation  is  equivalent  to  the  actual  situation;  for  if  we 
refer  it  to  coordinate  axes  moving  with  a velocity  -7  we  have: 

Figure  lid 

velocity 
P « V 
f 

stream  velocity 
at  oo  ■ o 


The  velocity  of  water  at  P is  still  v^,  the  computed  velocity.  Thus 
the  theorem  is  proved.  It  should  be  emphasized  that  the  velocity 
ccaaputed  is  steady,  but  that  the  velocity  observed  is  transient.  These 


two  velocities  can  be  ccmpared  according  to  the  above  theorem  only  when 
the  position  of  the  distribution  relative  to  the  point  P is  the  same 
as  the  position  of  the  ship  relative  to  P.  Since  the  ship  is  moving 
this  occurs  only  instantaneously. 
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Appendix  II 

Closed  Streamsurfaces  of  Point  Source  Distributions 


Rankine's  method  for  finding  the  limiting  streamlliies  is  not 
applicable  to  a problem  like  that  trtjated  here;  for  our  use  of  a 
double  source  and  a double  slnkc.  one  set  above  ^md  one  below  the 
water  surfaces  destroys  axial  symmetry,  and  a Stokes  stream  function 
does  not  exist.  Nor  Is  it  proper  to  neglect  one  set,  since  both  are 
equally  far  from  the  water  surface. 

Nevertheless,  It  Is  desirable  to  know  the  enclosing  streamlines, 
and  we  have  computed  them.  The  metnod  used  is  explained  in  this  Appendix. 

The  first  aim  is  to  find  the  curves  on  which  v^.  k o and  Vg  * o 
In  the  plane  of  the  sources.  The  intersections  of  these  curves  are 
either  sources  or  stagnation  points.  At  a stagnation  point  stream- 
lines meet  at  right  angles  For  a stagnation  point  in  froi:-t  of  a 

source,  one  of  the  two  Incoming  streamlines  is  from  the  source,  the 
other  from  the  stream.  The  two  outgoing  streamlines  form  the  longituuirMs^I 
midsection  of  the  Ranklne  solid  (see  Figs-  14,  20  and  21).  Also, 
streamlines  must  cross  the  Vx  » o curves  vertically  and  tne  Vg  “ o 
curves  horizontally.  Thus,  once  the  vx  and  Vg  « o contours  are  obtained 
the  streamlines  may  be  sketched  with  conaldei*able  accuracy. 

The  equation  which  the  Vj;  contours  satisfy  is  In  genera].  >\omogsneoi 
so  it  is  easier  to  solve  than  the  equation  for  v^^  ® o,  which  lai  not. 

For  our  problem,  the  v^^  » o curve  may  be  obtained  by  finding  the  locus 
of  points  for  which 


V, 


O) 


0) 


-\r  ' ■ _ . 2^- 


0 


f»  50  ■” 


( 


I 


> 


/ ci>  \ i ■>.  / r \ 

[ ;V  ; = - 12  V. 


We  thus  Ignore  the  sinks,  but  their  effect  Is  negligible  near  the  sources, 
For  an  arbitrary  planar  distribution  of  n sources,  the  contours 
of  Vy  satisfy  exactly  the  equation: 

(24) 


>y  v;. 

The  labor  of  solving  this  equation  can  be  greatly  reduced  by  consldez»~ 


ing  instead  the  n equations  of  the  form 

rlL4i.  — V - ' 


i = {25) 


where^  H . The  sun  of  these  eqs.  Is  £q.  (24)^  so  that  If  we 

solve  each  of  them  for  at  the  same  point  P and  add  ve  have  the 
total  v^^o  But  since  the  equations  are  all  of  the  same  form  ve  need 
solve  only  one  of  them,  say  the  J,th,  for  the  contours  of  Then, 

using  the  fact  that  for  points  Pj^  and  Pj  with  the  ssme  coordinates 
relative  to  sources  and  mj  respectively. 


(see  Figs  12a),  we  can  evaluate  v^  at  a particular  point  by  a simple 
process  of  multiplication  and  addition  (see  Flgo  12b). 


Fl{pire  12a 


Figure  12b 
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It  Is  alwayo  helpful  to  know  what  to  expect.  Stag;.iation  points 


on  a piano  of  symmetry  can  he  found  more  easily.  In  our  case.,  the 
equation  for  stagnation  points  on  the  X-axis  is: 

~7T-  - V =0 , 

and  it  was  found  that  there  are  2,  i,  or  0 stagnation  points  for  f 
less  than,  equal  to,  or  greater  than  soxae  critical  f . This  critical  f 
is  0.6204  ( )^.  For*  = 100  ft^  it  is  approximately  8.8  feet. 

The  streamlines  for  the  different  cases  ai*e  illustrated: 


Figure  13a 


Figure  13h 


Plg?jre  13c 


Since  f s 10  in  our  problem,  this  tells  us  our  streamlines  are  like  those 
in  Fig.  13c. 
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Appendix  III 


Pressure  Beneath  a Moving  Ship 
Thoughtless  application  of  Bernoulli -s  equation  leads  to  an 
instructive  paradox,  which  we  will  discuss  briefly.  T5ie  pressure  wl?vl 
be  calculated  by  two  methods,  with  apparently  different  results:  A. 
Consider  the  actual  situation,  depicted  in  Fig.  14.  The  ship  moves 
with  a velocity  V,  and  the 
water  far  fr<HQ  the  ship  is 
stationary.  At  the  points 
a and  b the  water  velocity 
Is  V. 


Bernoulli's  formula 
(with  the  obvious  p gh 
omitted)  asserts  that 


Figure  14 


<5Q 


i- 


Since 


0 — i 

r “x7 


y V 


\COl 


AT..  =r  a. 

’ V"  - V 

B.  Consider  now  the  fictitious  situation  In  which  the  ship  is  at  rest 
and  the  watet*  moveg^past  It  with  velocity  - V.  The  points  a and  b are 
aovj  stagnation  points,  and  the  water  velocity  at  is  '>V.  Indeed,  if 
we  denote  velocities  in  this  situation  by  v*  , in  the  former  by  v,  we  have 

- V 

Bernoulli's  formula  is 


COKFIDSmAL 
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cciTTiDEmm. 

- 33  ^ 


whence 

f Poo  + “ i ? 


(27) 


The  dlsorepaney  hetween  (26)  and  (2?)  is  apparent. 

To  resolve  it  we  note  that  the  nse  of  Bernoulli's  foxioala  in  A 
is  inproper,  since  the  pressure  at  a stationary  point  under  a moving 
ship  changes  in  timet  hence  the  situation  is  not  a stead;7  state. 

B does  represwit  a steady  condition,  but  not  the  actual  state  of  affairs, 
nevertheless,  it  leads  to  the  correct  pressure,  as  will  now  be  shown. 

For  a non*steady  state,  Bernoulli's  equation  must  be  replaced  by 
Buler's  idiich,  for  irrotational  flow  of  an  incoBqpressible  fluid,  takes 
the  form 

— l.,_const. 

J ^ 

where  ^ is  the  velocity  potential.  Mow 

<j>  » i^  .Z)  wKere  A — Vt  ; 


hence 


On  Inserting  this  in  Euler's  equation,  we  have 

(p  + rxV)^  = P '*’’1 

or,  since 

p * p-VxV  - ^ (28) 

which  is  Eq.  27* 

The  term  - 1/2 j?v^  makes  a relatively  small  contribution  to  p 
Hence  a knowledge  of  v^  is  generally  sufficient  to  understand  the 
pressure  field  of  a moving  vessel. 


Poo  • 
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In  the  absence  of  significant  flow  data,  we  can  test  our  theory 
only  by  computing  pressures  via  fonaula  (24),  which  reads  in  practical 
units 

A liAr.Wtoc  «X  V.V.W Wi  7:^ L ) 

129) 


Ap  (m  \vr-Vies  otvuaW)  = l.o<»\/  ^ 


The  ship's  velocity  (relative  to  the  water  at  ooi),  V,  is  hex*6 
expressed  in  knots,  and  Ap  increase  in  pressure  at  the  bottom 

over  its  normal  value,  j)gh.  ^cept  for  the  correction  term  propor- 
tional to  v2,ap  has  the  sign  of  v^.  This  agrees  with  obseinratlons . 

For  numerical  comparison  we  choose  seme  data  on  pressure  distrl*' 
bution  contained  in  NOL  Report  No.  504.  S.S.  Robin  Tuxford,  having  a 
length  of  48o  ft.  and  a beam  of  66  ft.  Is  closest  to  the  dimensions  of 
the  theoretical  model  among  the  four  boats  for  which  measurements  are 
listed.  It  developed  a maxisam  pressure  reduction  of  17  inches  in  the 
keel  section  under  the  fore  part  of  the  boat  when  proceeding  at  10,2 

knots,  and  the  depth  of  the  channel  was  40  ft.  under  these  conditions 

'Vi. 

Pig.  24  shows  tliat  * 0.135,  “ 0.  Hence  the  linear  term 

V » 

in  EQL.  (25)  contributes  AP  * " 14.9  the  other  - i.o  in.  The  total 
of  l6  in.  compares  well  x^lth  the  measured  value.  The  agreement  is  not 
Kiuch  worse  for  the  other  vessels  discussed  in  KOL  Report  Ho.  504,  even 


though  their  dimensions  differ  somewhat  from  those  of  our  model. 
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?a)  Page  2,  WoP,  V?hitmore,  NOL  Report  No.  504,  "Representation 
of  ships*  pressure  fields  by  source-sink  models":  14  pp.,  11  plates; 

7 April  1942. 

Cb)  Page  3*  For  references  to  papers  developing  and  applying  various 
mathematical  ship  models  see  the  bibliography  cm  pages  35*^37  of 
R.  Quilloton.  "Potential  Theory  of  Wave  Resistance  of  Ships  with 
!l^le8  for  its  Calculation",  Society  of  Naval  Architects  and  Karine 
Engineers,  Hew  York,  advance  paper  No.  2 for  the  6 and  7 September 
1931  meeting. 

(c)  Page  3.  See  pp.  407-413  of  L.M.  Hilne-Thomson,  Oeoretical 
Hydrodynamics . 2nd  ed.;  600  pp.,  1949*  This  book  is  the  primary 
hydrodyxmmiCB  reference  for  this  paper » 

(d)  Page  9.  This  equivalence  is  usually  left  Implicit.  A proof  is 
given  Appendix  I,  p.  27. 

(e)  Page  9.  Placing  the  source  and  sink  below  the  surface  rather  than 
in  it  was  suggested  by  J.K.  Limde,  "On  the  linearized  theory  of  wave 
resistance  of  displacement  ships  in  steady  and  accelerated  sction”, 
Society  of  Ka.va.1  Architects  and  Marine  Engineers,  New  York,  advance 
i>aper  No.  1 ror  the  6 and  7 September  195^  meeting.  This  nay  not  be 
the  best  placement.  Its  Interesting  effect  on  the  Rankine  solid  is 
discussed  in  Appendix  II,  p.  29*  Its  effect  on  bottom  velocity  cannot 
be  great  for  the  situation  we  analyze. 

I£i  Page  15.  The  method  by  which  these  streamlines  were  found  is 

presented  in  Appendix  II,  p.  29. 
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Cg)  Pa(56  21c  E,  Kadeluag.  "Dae  elektrleche  Peld  In  Systeraen  von 
regelmaue&sig  angeordneten  Punktladungen”  ^ Phyolkallsche  Zeltschrlft. 
12.  524-533  C1913)^ 

(h)  Page  2i,  H»Hargenau  and  OoM.  Murphy,  The  Mathematieg  of 
and  Chemistry.  1943.  P«  115. 

ti)  Page  29o  Milne-lhacpson,  op.cltc . pp. 110-111. 
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x(in  feet) 


vs.  X for  two  sources  m and 

Disrance  between  sources  354  ft. 
Distance  between  zeros  360  ft 
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